Bifurcation diagram and pattern formation in superconducting wires with electric currents 
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Abstract. We examine the behavior of a one-dimensional superconduct- 
ing wire exposed to an applied electric current. We use the time-dependent 
Ginzburg-Landau model to describe the system and retain temperature and ap- 
plied current as parameters. Through a combination of spectral analysis, asymp- 
totics and canonical numerical computation, we divide this two-dimensional pa- 
rameter space into a number of regions. In some of them only the normal state 
or a stationary state or an oscillatory state are stable, while in some of them two 
states are stable. One of the most interesting features of the analysis is the ev- 
ident collision of real eigenvalues of the associated PT-symmetric linearization, 
leading as it does to the emergence of complex elements of the spectrum. In 
particular this provides an explanation to the emergence of a stable oscillatory 
state. We show that part of the bifurcation diagram and many of the emerging 
patterns are directly controlled by this spectrum, while other patterns arise due 
to nonlinear interaction of the leading eigenf unctions. 

PACS numbers: 74.20.de 74.25.sv 74.25.dw 

We consider a one- dimensional superconducting wire of finite extent. An electric current 
is fed into the wire at one of its ends creating a voltage difference across the wire. This is 
a canonical problem that has received considerable attention since it involves the case of a 
resistive state in which a normal current and a superconducting current coexist. One of the 
intriguing phenomena associated with this state is the formation of phase slip centers (PSC). 
These are points in space-time where the order parameter in the time-dependent Ginzburg 
Landau equation (TDGL) vanishes. In fact, as was pointed out by Ivlev and Kopnin [T], phase 
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slip centers can be thought of as vortices in space-time. The appearance of phase shp centers 
is related to oscillations found numerically through the emergence of time-periodic solutions. 
The phase slip centers and the associated oscillations can be indirectly observed experimentally 
via the appearance of steps in I-V curves ([2], [3], Another type of behavior found in the 
resistive state involves stationary solutions of the TDGL [2]. In this case the gauge invariant 
quantities reach a steady state. 

One goal of this paper is to understand the origin of the different patterns observed in the 
resistive state. Another goal is to compute the bifurcation diagram in the parameter space. In 
particular we will explain why and when oscillatory solutions emerge. We will also consider 
the loss of stability of these oscillatory solutions as the temperature is lowered below a critical 
value that depends on the applied current /. The key idea is that the oscillations appear as 
a consequence of a Hopf bifurcation driven by a PT-symmetric spectral problem. A crucial 
role in the analysis is played by the dependence of this spectrum on the applied current. An 
additional goal is to elucidate the appearance of hysteresis in I-V curves in the present setup. 

Our starting point is the time-dependent Ginzburg Landau model that we write in a nondi- 
mensional form: 

ipt + i'-pi^ = i^xx + - \i^\^ip- (1) 

Here ip is the complex-valued order parameter, ip is the electric potential and F is proportional 
to Tc — T. Conservation of the current / implies the relation 

i^rx - i'xr) - <y^x = I, (2) 

where a models the Ohmic resistivity. (In equations ([T]), and all subsequent equations, we 
use a variable in subscript to denote a partial derivative.) The wire is assumed to extend along 
—L < X < L, and it is assumed that tp{±L,t) = 0. The main conclusions below are also valid 
for other boundary conditions such as ip^i^Lji) = 0. In order to concentrate on the main 
features of the phase transition mechanism, we take a = L = 1. This enables us to concentrate 
on the key parameters I and T. Some TDGL models include a factor + C2 1 P in front 

of the left hand side of equation ([1]). We deal here with the small (2 limit, but our essential 
results are valid for finite positive ^2- 
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Figure 1: The real parts of the first 6 eigenvalues of the PT- symmetric spectral problem 

To understand the phase transition from the normal state to the superconducting state we 
linearize the TDGL ([T]) about the normal state ip = 0, ip = —Ix. Writing ip{x,t) = u{x)e^^~'^^^ , 
we obtain for u{x) the spectral problem 



The spectral problem ([3]) is called PT-symmetric, since it is invariant under the joint transfor- 
mation of X —X (parity) and complex conjugacy (time reversal). The normal state thus loses 
its stability when F > Real(A(/)). However, since the spectral problem ([3]) is not self-adjoint, it 
is not clear at all that the spectrum is real. On the other hand the PT symmetry provides some 
useful information on the spectrum. Spectral PT-symmetric problems have attracted some 
interest in recent years following the numerical observation of Bender and Boettcher [6J that 
the spectrum of certain PT-symmetric problems is real. While ref. [0] considered a problem on 
the entire line, we deal here with a problem in a finite interval. When 1 = the spectrum is of 
course real. The PT-symmetry implies that if {X,u{x)) is a spectral pair with complex A, then 
also (A*,M*(— x)) is a spectral pair. Since the spectrum depends smoothly on / as long as the 
eigenvalues remain separated [7], a real eigenvalue cannot split spontaneously into a complex 
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pair. This implies that at least for small / all eigenvalues are real. However, when the current 
I is large, the lowest eigenvalues (in absolute value) are shown to satisfy A = 0{H), namely, 
to leading order they are purely imaginary. This implies that eigenvalues indeed collide as / 
increases. Specifically we find that the first such collision occurs when the first and second 
eigenvalues approach each other and collide at a critical value Ico ~ 12.31. 

At the collision point, the geometric multiplicity of the eigenvalue is 1. To find the behavior 
of the spectrum near Ico we set the current / to be / = Ico + eo,. Here e is a small positive 
number, and we introduce a to determine through its sign the direction in which we move from 
Ico- We then consider an expansion of the form 

A = /io + e^^^/Ui + + u = uq + e^^'^ui + eu2 + .... (4) 

The nonanalytic nature of the expansion for A is a consequence of the Jordan form of the 
spectral problem at the critical value / = Ico- The leading order term in (jlj) is found to be 
fiQ ~ 0.71, with an associated eigenfunction uq that we normalize by Uq{0) = 1. The first order 
correction /xi is conveniently expressed through the auxiliary function K{x) that solves 

K^^ + ixIcoK + i^oK = Mo, ^(±1) = 0. (5) 

Writing uq = Ur + iUi, and defining ai = 2 j\xUrUi dx and b = J\Kuo dx, one obtains 
fil = —aai/b. A numerical integration gives ai ^ 0.29 and b ^ 0.12. Since ai/b > 0, we see 
that when a < 0, i.e. when J is a little smaller than Ico, there are two real solutions; these 
are the first two real eigenvalues just before the collision. However, for / beyond Ico, that is 
for a > 0, the single eigenvalue /iq splits into a pair of complex eigenvalues. It can be further 
shown that fi2 is a single real number, i.e. it is the same for both splitting eigenvalues [S]. The 
analysis above shows that the real part of the leading eigenvalue is not an analytic function of 
the current at / = Ico- In fact, its derivative blows up as Ico is approached from below. On the 
other hand, the real part of the first eigenvalue (pair) is a smooth function of I just above Ico. 
This analysis holds for any later collision of real eigenvalues as well. It is in agreement with 
the numerical calculation presented in Figure [H 
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We computed the first few eigenvalues numerically as they increase past special collision 
points. Increasing I beyond Jco, the first two eigenvalues move as a complex pair according to 
the PT-symmetry. The real parts of the first six eigenvalues as a function of / are plotted in 
Figure [TJ We see there that respective pairs of eigenvalues collide at successive critical values 
of /. 

The normal state becomes unstable at that value of F for which F — Real(A) = 0. For 
/ < Ico the first eigenvalue A(/) is real. When the temperature is sufficiently low, i.e. when 
F = A(J), the normal state loses stability. Proceeding to high order terms in the bifurcation 
expansion it is found that the bifurcation branch that emerges at F = Fi(/) = A(/) is stable 
for I < Ik ^ 10.92. In this regime, i.e. when / < Ik and F > F(/), the bifurcating solution 
converges to a stationary solution. By 'stationary' here we mean that writing ip = /e*^, the 
gauge invariant quantities f{x,t), q{x,t) = Xx{x,t) and 9{x,t) = Xt{x,t) — (p{x,t) converge to 
stationary functions fo{x), qoix), 9o{x). Once / crosses the critical collision value Ico and the 
eigenvalue splits into a conjugate complex pair, the phase transition temperature is determined 
by the condition F = Real(A(/)) = Fi(/). Thus, for / > a Hopf bifurcation occurs and the 
solution to the full TDGL is periodic. 

Consider now a current / > Ico- When F is below Fi(/), the positive real part of the 
spectrum dominates, and the normal state is stable. Increasing F with I fixed we see that 
when F = Fi(/) a Hopf bifurcation into a periodic solution takes place as explained above. In 
addition to determining the bifurcation curve F = F(J), the spectral problem ([3]) can also be 
used to compute the bifurcating branch, which is always stable, in the periodic regime. To 
see this, fix a current / greater than the critical value Ico- Let the ground state of equation 
consist of the eigenvalue pair A^ ± ^A,, with associated eigenfunctions Wi{x) and W2{x) 
related by ^1(0;) = W2{—x). We normalize both eigenfunctions by the condition Wi{0) = 1. Set 
the temperature to be slightly below the critical value determined by F = Fi(/), by selecting 
F = A.r + e^. Neglecting a short time interval during which transients decay, the asymptotic 
solution of the full TDGL ([I])-(l2l) is found to be of the form 

u{x,t) = eA (exp {i{uje'^ + Xi)tj wi{x) + exp (-i(cje^ + A^) t) W2{x) + O(e^). (6) 
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The amplitude A and frequency uj are constants that are determined by the ground state Wi 
and W2- They are computed numerically for each current /. For instance, when / = 20 we 
found A = 0.921, uj = 1.8. One can draw a number of conclusions from the expression ([6]). 
First, the period of the oscillations is not exactly the imaginary component of the eigenvalue, 
but rather it is has a correction due to the nonlinear interaction of Wi and W2- Secondly, the 
solution at X = is ^(0, t) ~ 2eA cos ((Aj + uje'^)t). Therefore we obtain a phase slip center that 
is periodic in time at x = 0. Ivlev and Kopnin pLj made the nice observation that a PSC can 
be thought of as a vortex in space-time. In this sense, the solution structure given in equation 

indicates that the PSCs constitute a periodic placement of degree-one space-time vortices 
with period P = n / {Xi + ue"^). The curve ri(/) is depicted by the solid line in Figure [21 

So far we have concentrated on the smooth bifurcation of the normal state into a periodic 
state or into a stationary state. It turns out, though, that there are regions in the parameter 
plane where two metastable states coexist. The transition between them is nonsmooth, and 
therefore it is associated with hysteresis. We already pointed out above that for Ico > I > Ik 
the normal state bifurcates into an unstable branch. This hints that the phase transition there 
is nonsmooth. Indeed, we identified another curve in the phase plane, that we call r2(/), above 
which the stationary state is stable. The curve r2(/) is depicted as a dashed line in Figure El 

To understand the loss of stability of the periodic state, we recall that the Hopf bifurcation 
that led to it was triggered by the normal current contribution to the potential term in 
equation ([T]). Near the transition curve F = ri(/) the magnitude of the order parameter 
is still small, and essentially the entire current is normal. As the temperature is lowered (i.e. 
F increases), \ip\ grows and so does the supercurrent, implying via ^ that the normal current 
decreases. This effectively returns the system to the small I regime where the bifurcation to 
a steady state is favored. We thus obtain a third bifurcation curve F3(/) where the periodic 
state loses it stability. 

At this point we make reference to Figure [2] and consider the different regimes in the (/, F) 
plane. The solid curve provides the critical temperature F = Fi(/) along which the normal state 
loses its stability. A stable stationary state exists above the dashed line that represents a second 
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Figure 2: The phase diagram of the different stable states in the temperature-current plane. 
The parameter F is proportional to Tc — T. The curves Fi(/), F2(/), F3(/) are drawn with 
solid line, dashed line and dotted line, respectively. The meaning of the different curves and 
regimes is explained in detail in the text. 

curve F2(/). For I < Ik the normal state bifurcates into a stable stationary superconducting 
state. For / > Jco, on the other hand, the normal state bifurcates into a state that exhibits time- 
periodic oscillations. When / > Jco, and the temperature is further lowered (F is increased), 
the periodic state loses its stability at a third critical temperature F = F3(/) represented by 
the dotted line in the figure. The curves F2(/) and F3(/) intersect at Iq. For / > J^, the 
curves F2(/) and F3(J) coalesce. The frame on the left depicts the bifurcation curves over a 
large (/, F) area, while the frame on the right concentrates on the interesting area near the 
point (Jco, Fi(Jco)), where Fi(Jco) ~ 7.11. The parameter plane is partitioned into 5 domains. 
In domains 1,4 and 5 there is a single stable state - the normal state in region 1, a stationary 
state in region 4 and a periodic state in region 5. In region 2 there are two metastable states - 
normal and stationary, while in region 3 a stationary and a periodic state are both metastable. 

We proceed to draw two further conclusions related to the bifurcation diagram. It is useful 
to do so in the context of I-V curves. These curves are measured or computed for a fixed F, while 
the current I is raised or lowered adiabatically. When this process cuts through the metastable 
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Figure 3: The I-V curve for T = 6.3 and / increasing. Notice that, although the periodic state 
does not exist for such temperature, and thus there is no PSC here, the I-V curve does exhibit 
a jump discontinuity at I ^ 12.57 

regions 2 and 3 in Figure [2] a hysteresis is expected in the I-V curve. While such a hysteresis was 
predicted a long time ago, we point out that it is not always observed experimentally [1]. As 
can be seen in Figure [2], the metastable regions are quite small, and therefore it requires careful 
tuning to pass through them. Another comment relates to the formation of PSCs. These points 
in space-time where {ipl vanishes are often associated in the literature with jump discontinuities 
in the I-V curve. However, this identification works only in one way, and not all such jumps 
imply the presence of a PSC. For instance, we depict in Figure [3] the I-V curve for F = 6.3 
and / slowly increasing. For this F one never crosses an area in the parameter plane where the 
periodic state is stable, and therefore there is no PSC. Nonetheless, the I-V curve exhibits a 
clear discontinuity at about / ~ 12.57. The actual rule for lack of smoothness in I-V curves is 
that a jump discontinuity indicates a nonsmooth phase transition, while a discontinuity in the 
derivative indicates a continuous phase transition. 

To summarize, using a combination of asymptotic expansions, spectral analysis and canon- 
ical numerical computation applied to the time-dependent Ginzburg-Landau model, we have 
presented a full analysis of the behavior of a one-dimensional superconducting wire exposed to 
an applied electric current. In particular, retaining temperature and applied current as param- 
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eters, we have decomposed this two-dimensional parameter space into regions of stabihty of a 
normal, stationary and oscillatory state. The collision of real eigenvalues and the consequent 
emergence of complex spectrum in the associated linearized problem provides the explanation 
for the Hopf bifurcation leading to the appearance of the oscillatory state and the associated 
phase slip centers. From the theoretical standpoint, it also reveals a physically significant set- 
ting where PT-symmetry does not lead to reality of the spectrum, in contrast to its common 
role [6], [9]. The boundary of the basin of attraction of the normal state has been given precisely 
in terms of the real part of the leading eigenvalue in this linearized problem. The boundary 
between the basins of attraction of the oscillatory and stationary states has been calculated 
near the triple point using asymptotics, and has been computed numerically beyond this. In 
so doing, we have identified small regions in the parameter space where hysteresis should be 
anticipated. Finally the asymptotic structure of the periodic solution bifurcating off of the 
normal state has been developed for I above the first collision value Ico and for F just above 
the real part of the first eigenvalue. This expansion reveals the period of the oscillations and 
location of PSCs along the x-axis 
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